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Abstract 


We present a study of associated J/b+7 production through double Pomeron 
exchange at energies reached at the Fermilab Tevatron and CERN LHC based 
on the Ingelman-Schlein model for hard diffractive scattering and the factor¬ 
ization formalism of nonrelativistic QCD for quarkonia production. We find 
that this process (p+p —> p+p+J/ip + 7 +A) can be used to probe the gluon 
content of the Pomeron and study the nature of hard diffractive factorization 
breaking. 
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In high-energy strong interactions the Regge trajectory with a vacuum quantum number, 
the Pomeron, plays a particular and very important role in soft processes in hadron-hadron 
collisions JIJ. However, the nature of the Pomeron and its interaction with hadrons remain 
a mystery. Ingelman and Schlein 0 pointed out that hard diffractive scattering processes 
would give new and valuable insight about the nature of the Pomeron. It is assumed that 
the Pomeron, similar to the nucleon, is composed of partons, mainly of gluons, and the 
hard diffraction can be studied in a factorized way: First, the Pomeron is emitted from the 
diffractively scattered hadron, then one parton of the Pomeron takes part in hard subpro¬ 
cesses. Therefore, the partonic structure of the Pomeron could be studied experimentally. 
The Pomeron has a partonic structure was confirmed subsequently by various experiments 


However, there is an important problem, the factorization problem. The Ingelman and 
Schlein (IS) model for hard diffraction is based on the assumption of hard diffractive scat¬ 
tering factorization. Recently, a factorization theorem has been proven by Collins 0 for the 
lepton induced hard diffractive scattering processes, such as diffractive deep inelastic scat¬ 
tering (DDIS) and diffractive direct photoproduction of jets. In contrast, no factorization 
theorem has been established for hard diffraction in hadron-hadron collisions. At large \t\ 
(t is the square of the hadron’s four-momentum transfer), where perturbative QCD applies 
to the Pomeron, it has been proven that there is a leading twist contribution which breaks 
the factorization theorem for hard diffraction in hadron-hadron collision ||. This coher¬ 
ent hard diffraction was observed by the UA 8 Collaboration in diffractive jet production, 
in this experiment t is in the region —2 < t < — 1 GeV 2 ||. On the other hand, for the 
Pomeron at small \t\ , nonperturbative QCD dominates, it is unclear whether hard diffrac¬ 
tive factorization is valid or not in hadron-hadron collisions. In phenomenology, the large 
discrepancy between the theoretical prediction and the Tevatron date on the diffractive pro¬ 
duction of jets and weak bosons, at ai, signals a breakdown of hard diffractive factorization 
in hadron-hadron collisions [9j], but the nature of hard diffrative factorization breaking is 
still unclear. 

In a previous paper, we have studied associated J/i^ + 7 production in single diffractive 
(SD) scattering at the Fermilab Tevatron and CERN LHC energies based on the IS model 
for hard diffractive scattering and the nonrelativistic QCD (NRQCD) factorization scheme 
for quarkonia production [|TD[. In this BRIEF REPORT, we complete our study of associated 


J/ifj + 'y diffractive production by discuss associated J/V’ + q production at large Pt through 
double Pomeron exchange (DPE): 


p(Pp) + P(Pp) - P(Pp) + P(Pp) + J/^P) + l(k) + X. ( 1 ) 

This process is of special interesting because the large Pt J/i) produced is in the central 
rapidity region and is easy to be detected through its leptonic decay modes, and the J/^’s 
large Pt is balance by the associated high energy photon. We will see the measurement of 
this process at the Tevatron and LHC would shed light on the nature of the Pomeron and 
the hard diffractive factorization breaking. Furthermore, this process is also interesting to 
the study of heavy quarkonium production mechanism, and the Pt smearing effects. 

Within the NRQCD framework [[H],0] J is described in terms of Fock state decom¬ 
positions as 

I j/v) = 0(1) Iccpsi 1 ’]) + OMIofpfls} 
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+0(v 2 )|cc[ 1 S ( f ) ]c/) + 0(v 2 )\cc[ 3 S[ Lh) ]gg) 

+0(v 2 )\cc[ 3 Pj 1 ' 8 ' l }gg) + • ■ - , ( 2 ) 


where the cc pairs are indicated within the square brackets in spectroscopic notation. The 
pairs’ color states are indicate by singlet (1) or octet ( 8 ) superscripts. The color octet cc 
states can make a transition into a physical J/^j state by soft chromoelectric dipole (El) 
transition(s) or chromomagnitic dipole (Ml) transition(s) 

(cc)[“+‘Lf ] - J/i,. (3) 


These color-octet contributions are essential for cancelling the logarithmic infrared diver¬ 
gences which appear in the color-singlet model calculations of the production cross sections 
and annihilation decay rates for P-wave charmonia, and for solving the i^ 1 and direct J/%b 

IH3IJI4H 


‘surplus” problems at the Fcrmilab Tevatron 


The NRQCD factorization scheme 


13 iias been established to systematically separate high- and low- energy scale interactions. 
Furthermore, NRQCD power counting rules can be exploited to determine the dominant 
contributions to various quarkonium processes (l^|. For direct J/ij) production, the color- 
octet matrix elements, (0|(Pg' / ' i: ’[ 3 S'i]|0) , (0|(Pg^’[ 1 S'o]|0) and (01 (9g ^ [ 3 Pj] |0) are all scaling as 
So these color-octet contributions to J/ijj production must be included for consistency. 
On the partonic level, associated J/4> + 7 production are composed of the gluon fusion: 


S 7 
m c V c 


gi(pi) +g(p2) ^7 + (ccM 3 ^, 3 ^], 

9i(Pi) +g(P 2 ) -> 7 + (ccJpSjV Pj 8) ] ■ (4) 


The quark initiated subprocesses (qq channel) are strongly suppressed and will be neglected 
further. The color-singlet gluon-gluon fusion contribution to associated J/^ + 7 production 
is well known [ |T 6 || : 

^ (singlet) = D_£| M(g + g —* 7 + (oS)( 3 Si I, ])| 2 

(5) 


where s = (p\ + P 2 ) 2 , (0|(9 1 //,</ ’[ 3 S'i]|0) is the color-singlet matrix element which is related to 
the lepton decay width of J/rp in NRQCD 




r(j/^ —>> j + i )m 2 c 
(2/9)7 t e 2 c a 2 


( 6 ) 


where e c = 2/3. 

The average-squared amplitude of the subprocess g + g —■> 7 + (cc)[ 3 5i ] can be obtained 
from the average-squared amplitude of g + g —> 7 + (cc) [ 31 Q by taking into account the 
defferent color factor. The result is 


%\g + 9 -*■ 7 + (cc)[ 3 ^ 8) ] -> 7 + W] 

at 

= l6^T S|M(s + 9 ^ 7 + ( “ )[35il>1)|2 

x D- (0 |oWp Sl ]|0>. (7) 
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The average-squared amplitudes of the subprocesses g + g —■► 7 + (cc)^^] and g + g 
7 + (cc)[ 3 Pj 8) ] can be found in ||T7| , 


% [9 + 9 -> 7 + (cc)^^] ^7 + J/A] 

at 

1 E|M( g + g ^ 7 + (cc)[ 1 S < 8) ])| 2 


db 

dt 


16ns 2 

X ^(OloTYSollO), 

8 m c 

[s + 9 -* 7 + (ccJpP; 81 ] — 7 + JH>] 


£E|M( g + g ^ 7 + (cc)[ s P < 8 >])| 2 
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x (OlO^pflJIO), 

8 m c 

where the heavy quark spin symmetry 

(0\O j 8 ^[ 3 Pj]\0) = (2J + 1)(0|(9 8 ^[ 3 P 0 ]|0> 


( 8 ) 


(9) 


is exploited. 

Now we consider the Pt distribution of J/A produced in process (|T[). Assuming hard 
diffractive factorization, the differential cross section can be expressed as 

dcr = fip/p(t;i)fip/p(£,2)fg/ip(xi,Q 2 )fg/ip(x2,Q 2 ) 

^[singlet + octetjd^d^doqda^dt, ( 10 ) 

dt 

where £ 1 (^ 2 ) is the momentum fraction of the proton (antiproton) carried by the Pomeron, 
fg/w>(xi,Q 2 ) is the gluon distribution function of the Pomeron, and /p/ p (£) is the Pomeron 
flux factor integrated over t, 


/p/p (0 = /_ i fo/p{£,t)dt, 

fjp/ P (U) = ^ e~ Mt) F 2 (t ) 

= K^~ 2a ^F 2 (t), 


( 11 ) 


where the parameters are chosen as ||T 8 


K = 0.73 GeV 2 , a(t) = 1 + 0.115 + 0.26 (GeV" 2 ) t, 
F 2 (t) = e AM (valid at \t\ < 1 GeV 2 ). 


( 12 ) 


For numberical predictions, we use m c = 1.5 GeV, A 4 = 235 MeV, and set the factor- 
iztion scale and the renormalization scale both equal to the transverse mass of J/‘ 0, i.e., 
Q 2 = m t = (m 2 + Py), where Pt is the transverse momentum of J/0. For the color-octet 
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matrix elements (0|C?8^[ 3 Si]| 0), (0|(Pg^’[ 1 ,S'o]|0), and (0|(Pg^[ 3 Po]|0) we use the values de¬ 
termined by Beneke and Kramer |IIJ from fitting the direct J/ip production data at the 
Tevatron ]K| using GRV LO parton distribution functions p0| : 


( 0 |e> 8 /v, [ 3 Si]| 0 ) = 1.12 x 10 ~ 2 GeV 3 , 

(0|O 8 j/ ^5 0 ]|0) + ^|(0|C>8 /,/ ’[ 3 Po]| 0) = 3.90 x 10 ~ 2 GeV 3 . (13) 

m c 

The matrix elements (0|C?8 // ^[ 1 >S , o]|0) and (0|C>8^[ 3 -Po]|0) are not determined separately, so 
we present the two extreme values allowed by Eq.(|T3|) as 

1 S ' 0 saturated case :(0|C>8^ [ 1 S' 0 ]|0)= 3.90 x 10 ~ 2 GeV 3 

(0|0 8 J/,/ ’[ 3 Po]|0)=0, 

3 Pj saturated case /OlOg^’/Po] |0) = 1.11 x 10 - 2 m 2 GeV 3 

<0|C>^[ 1 S' 0 ]|0)= 0. (14) 

We use the hard gluon distribution function for the Pomeron ( 6 j: 

xfg/ip(x, Q 2 ) = fg 6 x(l - x), fg = 0.7. (15) 

We neglect any Q 2 evolution of the gluon density of the Pomeron at present stage. As usual, 
in order to suppress the Reggon contributions, we set £i ,£ 2 < 0.05. 

With all ingredients set as above, the Pt distribution of J/ip can be calculated in a 
standard way from Eq. (|10|). In Fig. 1 we show the Pt distribution Bda/dPx for associated 
J/^ + 7 production through DPE in pp collisions at a/s = 1.8TeV , integrated over the 
pseudo-rapidity region — 1 < p < 1 (the central region). Where B = 0.0594 is the J/ip —> 
p + p~ leptonic decay branching ratio. The lowest solid line is the color-singlet gluon-gluon 
fusion contribution, while the lowest dashed, dotted lines are 1 S'o 8 ' ) -saturated, and 3 S[ 8 ' > color- 
octet contributions respectively. The 3 Pj 8) -saturated color-octet contribution is almost the 
same as the 1 5 , g 8) -saturated case, so we don’t show it in the figure, and we don’t consider 
it further. For comparision, we also show the Pt distribution of J/ip produced in the 
inclusive process p + p —» J/-0 + 7 + X and the single diffractive process p + p—>p + 
J/i[} + 7 + X in the same kinematic region, the results are shown as the upper and middle 
lines of Fig. 1. The code for the lines is the same as the DPE case. As shown in Fig. 1, 
the color-octet 3 S [ 8 ' ) contribution is strongly suppressed compared with the others over the 
entire Pt region considered. The 1 S , q 8) -saturated contribution is smaller than the singlet 
contribution where Pt < 5 GeV. Integrated over the Pt region (2.0 < Pt < 5.0GeV), the 
color-singlet, 1 5'o 8 ) -saturated, and 3 Sj 8 ^ contributions to Ba DPE are 0.47, 0.12, and 1.4 x 10 ~ 2 
pb respectively. The total cross section times the leptonic decay branching ratio of J/pjB 
at a /s = 1.8TeV, integarated over the same Pt and r/ region for DPE, SD, and inclusive 
production of associated J/ip + 7 are 0.60, 19, and 4.0 x 10 2 pb respectively. The ratios of 
the total DPE and SD cross sections to that of inclusive production in the central region 
are i? DPE = 1.5 x 10 -3 and R SD = 4.8 x 10 ~ 2 respectively. 

I 11 Fig. 2 , we show the P T distribution of J/^>, Bda/dP T , integrated over the same 
pseudo-rapidity region at the CERN LHC energy y/s = 14 TeV for the inclusive, SD, and 
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DPE processes, the code for the lines is the same as Fig. 1. Integrated over the Pt region 
(2.0 < P T < 5.0GeV), the color-singlet, 1 S'®-saturated, and 3 Sj 8) contributions to Ba DPE 
are 1.2, 0.30, and 3.6 x 1CT 2 pb respectively. The total cross section times the leptonic 
decay branching ratio of J/^B at y/s = 14TeV, integarated over the same Pt and rj region 
for DPE, SD, and inclusive production of associated J/ip + 7 are 1.6, 88, and 3.1 x 10 3 
pb respectively. The ratios of the total DPE and SD cross sections to that of inclusive 
production in the central region are i? DPE = 5.2 x ICE 4 and R SD = 2.8 x ICE 2 respectively. 
We have varied the color-octet matrix elements (01[ 1 5' 0 ] 10) + ^(0\O^[ 3 P 0 ] |0) and 

TYl c 

(O|0 8 J/0 [ 3 Si]|O> by multiplied them by a factor between 1/10 and 2, the ratios J? DPE and 
R SD above are unvaried. This character demonstrates that the ratios f? DPE and R SD are 
insensitive to the values of color-octet matrix elements. 

In TABLE I., we show the ratios 


R dpe (P t ) 


R sd (P t ) 


da DPE / d(7 Ir,clusive 
dP T ' dP T 


da SP> j dcr Inclusive 

dP T ' dP T 


(16) 


at yfs = 1.8 and 14 TeV in the central rapidity region. 

From this table, we can see that R sp> (Pt ) is almost constant for the Pt region considered, 
so R SP> is insensitive to the Pt smearing effects. But R dpe (Pt) varies rapidly with Pt, from 2 
GeV to 5 GeV, which decreases by almost one order, so R DPE is sensitive to the Pt smearing 
effects, this can be attributed to the limited production phase space and an additional hard 
gluon distribution of the Pomeron in the DPE compared with the SD case. 

I 11 the above calculation, we use the standard Pomeron flux factor ( valid at small t ) that 
appears in the triple-Pomeron amplitude for SD. In order to preserve the shapes of the M 2 
and t distribution in soft single diffraction and predict correctly the experimentally observed 
SD cross section at all energies in p — p collisions, Goulianos 
the Pomeron flux in an energy-dependent way: 


proposed to renormalize 




(17) 


the renormalization factor D is defined as 


with 


D = min(l, —) 


/*£ max p 0 

N= dd dtfjp /p (t,t), 

^£min J—OO 


(18) 


(19) 


where £ m in = M^/s with Mg = 1.5GeV 2 (effective threshold) and £ max = 0.1(coherence 
limit). For yfs = 1.8,14 TeV, the renormalization factor D = 0.11,5.2 x 10 -2 respectively. 
Furthermore, by comparing the CDF data on diffractive W production with predictions 
using the Pomeron structure function measured at HERA, Goulianos concluded that the 
breakdown of hard diffractive factorization in hadron-hadron collisions is due to the break¬ 
down of the Regge factorization already observed in soft diffraction [p]]. In this picture, the 
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renormalization factor D indicates the factorization broken effects and takes the role of the 
survival probability for hadron emerges from the diffractive collision intact [^]. Using the 
renormalized Pomeron flux factor, the Px distribution calculated above must be multiplied 
by a factor D 2 and D for DPE and SD case respectively, therefore, the ratio i? DPE and R SD 
are proportional to D 2 f 2 and Df g respectively , hence they are sensitive to the gluon frac¬ 
tion of the Pomeron f g and the renormalization factor D which indicates the factorization 
broken effects. From other diffractive production experiments, f g can be determined, the 
renormalization factor D can be determined precisely from those ratios, and vice versa. So 
measuring these ratios can probe the gluon density in the Pomeron and shed light on the 
nature of hard diffractive factorization breaking. Furthermore, sence i? DPE is sensitive to 
the Px smearing effects, the experimental study of associated J/A + y production through 
DPE can give valuble information about the Px smearing effects which are needed to solve 
the large-z discrepancy seen by comparing NRQCD predictions with the HERA data on 
inelastic J/ip photoproduction i|?o| . 

Experimentally, the nondiffractive background to the diffractive associated J/ip + y pro¬ 
duction must be dropped out in order to obtain useful information, this can be attained by 
performing the rapid gap analysis or using the Forward Proton Detector. 

In conclusion, in this BRIEF REPORT we have shown that diffractive associated J/k’ + y 
production at large Px is sensitive to the gluon content of the Pomeron and the factorization 
broken effects in hard diffraction. Although the diffractive and inclusive production cross 
sections are sensitive to the values of coler-octet matrix elements, the ratioa of the diffractive 
to inclusive J/^+y production are not so and proportional to D 2 f 2 and Df g for DPE and SD 
case respectively, hence they are sensitive to the gluon fraction of the Pomeron and the the 
factorization broken effects. So experimental measurement of these ratios at the Tevatron 
and LHC can shed light on the nature of Pomeron and hard diffractive factorization breaking. 
Furthermore, sence the ratio of the total DPE cross section to that of inclusive production is 
sensitive to the Px smearing effects, the experimental study of associated J/^ + y production 
through DPE can give valuble information about the Px smearing effects. 
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Figure Captions 


Fig. 1. Transverse momentum of distribution Bda/dPx , integrated over the 

J/i/) pseudo-rapidity range rj < 1 (central region), for associated J/ijj + 7 production 
through double Pomeron exchange (lower), in single diffractive (middle), and inclusive 
(upper) processes at the Tevatron (yfs = 1.8 TeV). Here B is the branching ratio of 
J/^ —> n + n~(B = 0.0594). The solid line is the color-singlet gluon-gluon fusion contribu¬ 
tion, the dashed line represents 1 Sq ’^-saturated color-octet contribution, the dotted one is 
color-octet contribution. The inclusive and single production data are multiplied by 
a factor 100 and 10 respectively. Fig. 2. Transverse momentum of J/^^Pt) distribution 
Bdcr/dP T , integrated over the J/tjj pseudo-rapidity range \r}\ < 1 (central region), for asso¬ 
ciated J/ip + 7 production through double Pomeron exchange (lower), in single diffractive 
(middle), and inclusive (upper) processes at the CERN LHC (v/i = 14 TeV). Here B is the 
branching ratio of J/r/? —> h + /j,~(B = 0.0594). The solid line is the color-singlet gluon-gluon 
fusion contribution, the dashed line represents 1 S'q 8 ^-saturated color-octet contribution, the 
dotted one is 3 S{ 8 ^ color-octet contribution. The inclusive and single production data are 
multiplied by a factor 100 and 10 respectively. 



TABLES 


TABLE I. The ratio R se> (Pt), R e>pe (Pt) as a function of Pp at the Tevatron and LHC energies 
in central region (—1 < rj < 1) using the standard Pomeron flux (D = 1.0). 


P T (GeV) 

2.0 

2.5 

3.0 

3.5 

4.0 

4.5 

5.0 

R sd (P t ) (V s = L8 TeV) ( 10_2 ) 

4.9 

4.9 

4.8 

4.8 

4.7 

4.7 

4.7 

R sd {P t ) (Vs = 14 TeV) (10“ 2 ) 

2.9 

2.9 

2.9 

2.9 

2.9 

2.9 

2.9 

R dpe (P t ) (y/s = 1.8 TeV) (10" 3 ) 

2.0 

1.7 

1.3 

1.0 

0.71 

0.46 

0.25 

R dfe (P t ) (\/i = 14 TeV) (10“ 4 ) 

6.9 

5.8 

4.5 

3.3 

2.2 

1.1 

0.23 
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